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Interplay between electron spin and orbital pseudospin in double quantum dots 

Sooa Park * and S.-R. Eric Yangt 
Dept. of Physics, Korea University 
(Dated: February 2, 2008) 

We investigate theoretically spin and orbital pseudospin magnetic properties of a molecular orbital 
in parabolic and elliptic double quantum dots (DQDs). In our many body calculation we include 
intra- and inter-dot electron-electron interactions, in addition to the intradot exchange interaction 
of 'p' orbitals. We find for parabolic DQDs that, except for the half or completely filled molecular 
orbital, spins in different dots are ferromagnetically coupled while orbital pseudospins are antifer- 
romagnetically coupled. For elliptic DQDs spins and pseudospins are either ferromagnetically or 
antiferromagnetically coupled, depending on the number of electrons in the molecular orbital. We 
have determined orbital pseudospin quantum numbers for the groundstates of elliptic DQDs. An 
experiment is suggested to test the interplay between orbital pseudospin and spin magnetism. 

PACS numbers: 73.63.Kv, 73.23.Hk 
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I. INTRODUCTION 

Nanometer-scale quantum dots have potential for 
many technological applications 1]. A useful model of 
these systems is a parabolic quantum dot. In this model 
two-dimensional electrons are confined in a parabolic po- 
tential. The eigenenergies of a parabolic quantum dot 
are given by 



E n = Hu>(n + 1), rt = 0, 1, • 



(1) 



where lu is the characteristic frequency of the parabolic 
potential. This dot is rotationally invariant about the 
symmetry axis of the parabolic potential. (Hereafter this 
axis will be called the z-axis). The energy states for n = 
0, 1, 2, • • • are labeled as 's', 'p', 'd' ... orbitals, as in the 
real atoms. Each of these states is (n+l)-fold degenerate. 
The 's' orbital is non degenerate with z-component of 
the angular momentum quantum number a — while 
the 'p' orbital for each dot is doubly degenerate with 
a = ±1. Parabolic dots form basis for understanding 
more realistic dots, for example, elliptic dots [3]. 

When two of these dots are coupled an artificial di- 
atomic molecule can be formedQ. These DQDs have 
potential application in quantum computing since they 
may provide a basic building block of solid state realiza- 
tion of a quantum computer 0, IE IE IE IE El- Drawing 
the analogy with the transition metal compounds[lfJ, 
we expect to find interesting physics in the molecular 
orbital originating from the degenerate levels of a quan- 
tum dot, such as the p state. The orbital degeneracy 
of a shell makes the physics rich in the transition metal 
compounds. To illustrate the physics we consider two 
sites each with a two-fold degenerate orbital. The de- 
generate orbitals are denoted by a and b. These states 
can be labeled with a pseudospin value: the state a with 
pseudospin 1/2 and the state b with pseudospin -1/2. In 



such a model, there are on-site repulsion energies U a , Ub 
and U a b, intra-atomic exchange energy J, and the tun- 
neling energies t a and tb- When two atoms each con- 
tains an electron and under the following special case of 
t a = tb = t, U a = Ub = U and U a b = U — J, the effective 
model Hamiltonian for t <C t7 can be written lO] as 



H, 



At 2 



eff 



U-2J 



— + S 1 • S 2 



1 - 

4~ T1 



T 2 



"w U ~ 1 



T 2 



(2) 
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Here, spin and pseudospin of the electrons on the j-th 
atom are denoted by Sj and r respectively. This 
particular Hamiltonian is rotationally invariant both in 
spin and pseudospin spaces, i.e., it has the symmetry of 
SU{2) x SU{2) 10]. In the groundstate of this Hamilto- 
nian spins in different atoms are ferromagnetically cou- 
pled while their pseudospins are antiferromagnetically 
coupled. The ground-state energy is —At 2 /(U— 2 J). This 
state has the total spin and pseudo spin quantum num- 
bers S = 1 and r = 0, and is 3-fold degenerate. The first 
excited state with energy —At 2 /U has S = and r = 1, 
and is also 3-fold degenerate. The second excited state 
has zero energy with S' = r = 0or5 = T = l (10-fold 
degenerate) . 

Magnetic and excited state properties of DQDs are also 
of experimental interest. Recently an experimental inves- 
tigation of orbital magnetism of DQDs containing about 
50 electrons was carried out^^- The single-electron tun- 
neling spectroscopy in a finite source drain volt age|l3| al- 
lows one to explore the excited states of DQDs 14]. This 
technique is used to detect the singlet-triplet transition 
in DQDs^E], which plays an important role in two-qubit 
quantum gates of quantum computing pj]. 

The physics of the molecular orbital made of the dou- 
bly degenerate 'p' states of DQDs is expected to be more 
complicated than that of the effective Hamiltonian de- 
scribed in Eq. (J2J. In DQDs, 'p' states can be labeled 
with a pseudospin value: the state with a plus value of the 
z-component of the angular momentum has pseudospin 
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1/2 and the state with a minus value has pseudospin - 
1/2. An electron has thus a pseudospin 1/2 on top of 
ordinary spin 1/2. Unlike transition metal compounds, 
since U a b — U rotational symmetry of f is broken, i.e., 
the SU{2) x SU(2) symmetry is broken in DQDs. How- 
ever the system is invariant under pseudospin rotations 
about the z-axis. It will be interesting to investigate 
how pseudospins and real spins of two dots couple in the 
molecular orbital made of doubly degenerate 'p' states 
of DQDs. The other issue is the orbital magnetism of 
a DQD0,0,^}- The orbital magnetism is related to 
the z-component of total angular momentum. When a 
DQD contains two electrons it is zero. When the total 
number of electrons in the molecular orbital N p is differ- 
ent from two the z-component of the total angular mo- 
mentum can take non-zero values and therefore orbital 
magnetism may take a finite value. For elliptic quantum 
dots, the rotational symmetry about the z-axis is broken. 



We assume the energy level spacing of the harmonic 
potential is much larger than the strength of the Coulomb 
interaction: ftiu >> U . Only the 'p' orbitals are thus rel- 
evant (the molecular orbital originating from the 's' or- 
bitals are fully occupied with four electrons and can be ig- 
nored). This can be satisfied in self-assembled DQDsjlJ 
where the level spacing is hcu ~ 40meV. We have ob- 
tained wavefunctions, energy splittings, and degeneracies 
of states both for parabolic and elliptic DQDs. In our 
many body calculation we include intra- and inter-dot 
electron-electron interactions, in addition to the intra- 
dot exchange interaction of 'p' orbitals. Since pseudospin 
rotational symmetry is broken not all the eigenstates of 
the Hamiltonian can be also eigenstates of the total pseu- 
dospin operator r 2 . However, we find that the ground- 
states of both parabolic and elliptic DQDs happen to be 
eigenstates of r 2 . We find for parabolic DQDs that, ex- 
cept for the half or completely filled molecular orbital, the 
groundstate spins in different dots are ferromagnetically 
coupled while orbital pseudospins are antiferromagneti- 
cally coupled: the total spin of the electrons is maximized 
while the total pseudospin is minimized for all numbers of 
electrons except for N p = 4, 8. For elliptic DQDs ground- 
state spins and pseudospins are either ferromagnetically 
or antiferromagnetically coupled, depending on the num- 
ber of electrons in the molecular orbital: spin ferromag- 
netism occurs for N p — 3 and N p — 5 while pseudospin 
ferromagnetism occurs for N p = 2 and N p = 6. An exper- 
iment is suggested to test the interplay between orbital 
pseudospin and spin magnetism. 



This paper is organized as follows. In Sec. II we de- 
scribe our model Hamiltonian. In Sec. Ill we discuss sym- 
metry properties of this Hamiltonian. The proper basis 
states of the unperturbed Hamiltonian are given in Sec. 
IV. Low-lying states of parabolic and elliptic DQDs are 
given in Sec. V and VI, respectively. Conclusions are 
given in Sec. VII. 



II. MODEL HAMILTONIAN FOR DOUBLE 
QUANTUM DOTS 

We describe our model Hamiltonian for double quan- 
tum dots. If the interdot exchange energy is ignored, 
Coulomb interactions for the parabolic dots can be de- 
composed as follows: 

U a , a , = [ /d 2 xd 2 ^Va(x)| 2 |^(x')| 2 «(x-x') 



V ataf = I I d 2 xd 2 xVaW| 2 ^a'(x' + d)| 2 V (x-x'-d) 



J 1 ,_l= / / d 2 a;d 2 a;Vl(x)^_i(x)^I 1 (x , )v5i(x') 

xu(x-x'), (3) 

where (p a is the wavefunction of an electron with the 
z-component of the angular momentum a and v is the 
Coulomb potential. U a ^ a i is the on-dot direct repulsion 
energy between two electrons, one with the z-component 
of the angular momentum a and the other with a! . V ata ' 
is the same as U a , a ' except that it is the interdot direct 
interaction. J\ t —\ = J is the orbital exchange energy be- 
tween two electrons in different orbitals of the same dot. 
Since the direct interactions are independent of a and a' , 
i.e., J7 a ,a' = U and V a ,a.' = V, we have the Hamiltonian 
for parabolic DQDs with no hopping between dots 



Ho 



U 

\ jot jaal 
+ C ]lcr C j,-l,(T' C ilo"' C j,-l,o- 

+V ni aa n 2 a>a> + Ei n jaa 



(4) 



where 



and the fermion operator Cj a 



( c ]acr) annihilates (creates) the electron with the spin a/2 
in the state of the z-component of the angular momentum 
a on the j-th dot. The summations are performed for 
j = 1, 2, a = ±1 and a = ±1 (f, 1). 

The hopping Hamiltonian between double dots is given 

by 



Fhop 



E (4 



,C2aa + h.C. 



(5) 



where t denotes the tunneling energy of an electron be- 
tween the opposite dots. It is assumed that the hoping 
of an electron between dots changes neither the angular 
momentum nor the spin. Furthermore, we assume that 
the hopping energy t is not dependent of the angular 
momentum of the electron. 

The confinement potential of each quantum dot can be 
assumed to be parabolic with a distortion. The distor- 
tion an electron experiences has an elliptic form and is 
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represented by 



III. SYMMETRY OPERATORS 



Voiiiptic(r) = 2 [ — — ) exy, 



(0) 



where e is a positive parameter describing the deviation 
from parabolicity. The distortion breaks the U(l) sym- 
metry and lifts the degeneracy of a parabolic dot. It splits 
the 'p' orbitals of E\ — 2fnu with the angular momenta 
a = ±1 into two states 

|-1>' = -L(|l>-i|-l>) ,hvE[ = (2tkJ-e) , 

|1 >' = -= (|1 >+i\ -1 >) ,for^ = (2fiw + e) , 
v2 

(7) 

where |a > denotes the state with the angular momen- 
tum a = ±1 and E[ is the single particle energy of the 
elliptic dot. The total Hamiltonian for the DQD with the 
elliptic potential can be written as 



H + Km 



ptlCi 



(8) 



where Ti = Ho + Vh op is the total Hamiltonian for the 
DQD with a parabolic potential and Vcinptic represents 
the Hamiltonian of the elliptic potential. The elliptic 
potential term in the second quantized form is 



Vkiiptic = (2fiw - e)?7_i + {2hw + e)rji 



(9) 



where rjp is the number operator of particles in the states 
\P = ±1 >'. It can be written as 



\ I E 4^ C 3^ + l j ^l ,U i 
\ jaa jaa 



.(10) 



In terms of pseudospin operators the elliptic potential 
can be written as 



^elliptic = 2er y + EiN p , 



where 



E 

3<r 



(11) 



(12) 



In general, the intradot interaction is greater than the 
interdot interaction since two dots are separated with a 
distance d each other. One can assume that U > J > 
V t. We consider elliptic DQDs that are significantly 
distorted from parabolic DQDs, so we assume that e t. 
Both the elliptic potential and the hopping terms will 
be treated as perturbations. The number of 'p'-orbital 
electrons N p that can be injected in the system is from 
to 8. For these numbers, the low energy states will be 
obtained in Sec. V for parabolic dots and in Sec. VI for 
elliptic dots. Before these, we will explain the symmetry 
operators of our Hamiltonian and the basis states for H.q. 



Orbital pseudospin operators Let f be the total orbital 

(13) 



pseudospin operator defined by 
1 



T = o( T + T - 



with 



ir v = E t i+ = E 



T x ITy 



E r ^- = E Eii/^ 

3 3 V °~ / 



3 L 



( C tlcr C 3~ la ' C %-\,o C 3--^,°) 



(14) 



These operators have the following properties. When a 
single dot contains one electron then Tj + \ 4J->j= I ft>j 
and Tj_| JJ->j= 0, where | f|->j and | -l|>j denote the 
states of the j-th. dot with the up pseudospin and the 
down pseudospin of the electron, respectively. If a sin- 
gle dot contains two electrons with opposite spins then 
Tj-\ 1NJ->j= | JJ4>j, where | M>i= | t>j ®| §>j- A l so 
when two electrons with the same spin occupy a single 
dot then the operation with Tj + or Tj_ gives zero due 
to the Pauli exclusion principle. Lastly in the case each 
dot contains one electron we find t + | ltJJ->= | iytf> and 
T-l M>= I H>, where | 1*NJ->= | ft>i ®| J|>2- 

In general, neither parabolic dot nor elliptic dot has the 
orbital pseudospin number t as a good quantum number 
since [Tio,T 2 ] ^ holds. This can be understood as fol- 
lows. Let us define the total orbital pseudospin of j-th 
dot Tj as 7? — {Tj+Tj- + Tj-Tj + )/2 + r? . The ferro- 
magnetic exchange term of TIq in Eq. Q can be written 
as Tif.ox. = - J\[N P - J2j( T 3+ T j~ + T 3- T 3+)]- I* can b e 
shown that this term and t 2 do not commute yielding 



[H f . cx .,r 2 ] = 2J^(-1)^ 1 [t j+ (t 1 ,-t 2z )t j+1 ,_ 

-Tj-{T1 Z - T 2z )Tj +li+ ] ^ (15) 

where t^,± = Tx t ±. 

We will show nonetheless that for some states 



[Wf.ex., t 2 



>=0. 



(16) 



If | ip > is a one electron state of j-th dot then this equa- 
tion holds since Ti.f. cx .\'^ >— and Hf.oxT 2 !^ >= 
(T 2 \ip > is also a one electron state). Also if \tp > has two 
electrons in a singlet orbital-pseudospin state of the j-th 
dot \Tj = 0,r jz = >= ^(| m>j -| U>j) then the 
commutator operating to this eigenstate vanishes. It is 
because the commutator in Eq. i|15[l has Tj + , Tj- or Tj Z in 
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every term and those numbers for the singlet state are all 
zero as follows: r j+ ^(\ U>j -\ M>j) = ^(1 U>j 
-\ n>j) = 0, tj-^Q U>j -I M>j) = ^(1 U>j 
-| U> 3 ) = and r jz ±(\ U>j -I U>j) = 0. Hence, 
if any dot of double dots has a singlet orbital-pseudospin 
states with either t\ = or t 2 = 0, the number of the 
total orbital pseudospin r becomes a good quantum num- 
ber. 

Meanwhile, the z component of the total orbital pseu- 
dospin commutes for parabolic dots, [H, t z ] — 0, so 
that r 2 becomes a good quantum number. For elliptic 
dots, however, it is not a good quantum number since 

[Elliptic, Tz\ -f- 0* 

Spin operators Let S be the total spin operator. Both 
parabolic and elliptic dots have spin rotational symmetry, 
i.e., [Htot,S] — 0. This symmetry can be verified using 
the spin operators 

S+ = S x + iS y = c ]al C ja,~l 
jot 

S- = S x — iSy — ^ Cja ^Cjai 
ja 

ja 

Dot pseudospin operators One can consider another 
pseudospin depending on which dot an electron belongs 
to[l8l ITgj . If the electron resides on the lst(2nd) dot, 
the dot pseudospin of it is given by 4(— |). With the 
similar definitions to the spin and the orbital pseudospin 
operators, we have 

n 2 = - + n_n + ) + n 2 z (is) 

with 

f2 + = tt x + iQy = c laa C 2aa 
a a 

ota 

= i i^laa^ota - 4 aa C 2a a) ■ (19) 
aa 

In general, the number f2 is not a good quantum num- 
ber, i.e. [Ho,ri 2 ] ^ since the term Hf. ex . does not 
commute with Q, 2 . From the operational calculations, 
however, we find the number Q is well-defined for some 
eigenstates. If \ip > is a one-electron state of the j-th dot, 
then H f . ex |^ >= 0. And H f . cx ft 2 |V> >= since fl 2 \ip > 
is also a single-electron state. Hence, 

[H f . cx ,O 2 ]|V>=0. (20) 

We also see that [Ho,Q z ] = but [Vh op ,^ z ] ^ 0, and 
hence £l z does not commute with the total Hamiltonian. 

Particle-hole transformation Both parabolic and ellip- 
tic dots have particle-hole symmetry. For the number of 



electrons in 'p' orbital N p > 4, one can straightforwardly 
obtain the energy spectra and the corresponding states 
from the results for the systems of N p < 4 via a particle- 
hole transformation without detailed calculations. Here 
we will prove this for the parabolic case. Let us consider 
first the vacuum state in the hole space. The maximum 
number of electrons that we can insert on the double dots 
is 8. In this case there exists one possible state, i.e. all 
orbitals are filled with particles. This fully occupied state 
is defined as a hole vacuum \0 >h which is related to the 
particle vacuum |0 > as 

|0 >h— c\ll c \l ,-l c l.-l,l c l.-l -l c 211 c 21.-l c 2,-l.l c 2,-l.-l|0 > • 

(21) 

The hole operators are defined as hj aa — Cj a<7 (hj aa = 

Cj aa ). They, of course, fulfill the anticommutation rela- 
tions 

\hjcxo~ 3 ^j'ot' ty 1 } — ^jj' ^aa'^aa' (22) 
\hj acr , hj'a'tx' J = 3 (23) 

and operating the hole annihilation operators on |0 >h 
gives zero, i.e. hj arT \0 >h= 0. The total number of holes 
is given by 

N h = J2 h ]ot,h jaa = 8-N p . (24) 

jotcr 

By replacing c jaa {c] arT ) with h] a(T {h ]aa ) and rear- 
ranging the operators, our Hamiltonian H can be easily 
transformed as 

H(U,J,V,t;N p ) = 

H h (U, J, V, -t\ N h ) + &Ei + (3C7 — J + 4V)(4 - N h ) . 

(25) 

In Eq. I|25|) . the hole Hamiltonian H h is defined by 

H h = H\ + < op (26) 

with 

\ ja jaal J 

+ J^2 h ]la h ]-l,a> h jl<?' h :i-l,°- 
jaa 

+V £ n h laa n h 2a , a , + E^J2 n 1 a , (27) 

olgol' <t' jaa 

where n h jaa = h' jaa hj a(T and E\ = -E t , and 

V t P = ~t E ( h \ a « h ^ + h - c -) • (28) 

The energy eigenvalues for 7i h with Nh holes is nothing 
more than the spectra for H. with the same number of 
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(a) 



B'(l) 


o=l 


ce=l 





FIG. 1: (a)Basis set B 1 (N P = 1). T (I) is the electron with 
up (down) spin. Upper (lower) half circle represents a state 
in the 'p' orbital with z-component of angular momentum 
a — 1( — 1). Not all the basis states are shown. Other states 
can be obtained by flipping all spins, by altering the signs 
of all angular momenta, or by doing both. (b)One of the 
eigenstates for parabolic dots is shown as an example. This 
state has the quantum numbers a = 1 and a — 1 and is given 
by a linear combination of the basis staes in the set B 1 (l). 



b'(2) 0=2 0=0 



11=2 



a-0 




FIG. 2: Some basis states in the set B 1 (N P — 2). States with 
two electrons in a single dot are not included in this basis set 
since we are interested only in the low energy states. 

electrons, which makes it possible that the low-lying en- 
ergies for N p > 4 are obtained from the results for N p < 4 
by using — t instead of t and 2V/, = 8 — N p . One can also 
have the low- lying states by retransforming hj aa , hj aa 

and |0 >h into c JQ(J , ct a and |0 >. The proof for the 
elliptic case is the same as for the parabolic one except 
that 

Wtat (U, J, V, t, e; N p ) = H h tot {U, J, V, -t, e; N h ) 

+16E 1 + (3U-J + W){4-N h ) (29) 

where Ht t = H h + V* niptic and 

KlUptic = " X! (/' - I - // '- 1 ' T h l i^'' 1 ") 

+E*N h . (30) 

Hence, all eigenstates and energies for elliptic dots for 
N v > 4 can be obtained from those for N v < 4. 



IV. BASIS STATES AND DEGENERATE 
PERTURB ATIONAL CALCULATION 

The total Hamiltonian of our system is Hq + V\ lop + 
Elliptic- In the parabolic DQD Veiuptic is absent and we 
will treat Vhop in degenerate perturbation theory using 
the eigenstates of Hq as basis states. In the elliptic DQD 
we will use the eigenstates of Hq to find eigenstates of 



(a) 

B'(3) a=3 o=l 




(b) 



B 2 (3) 




a=3 




ee-ee 


a=l 




ee.e© 



(c) 



B 3 (3) 


0=1 


0=1 


*<®e-©0>- *<e©-e©> 



FIG. 3: Some basis states in the three sets from B 1 (N P = 3) 
to B 3 (N P = 3). 
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B'(4) 
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(b) 



B ! (4) 




o=0 




e©.@© 





B 3 (4) 




o=0 






S3 


a=0 


*<©©-©©> 


©©■ ©© 
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(d) 



B 4 (4) 


o=0 


a=2 


*<©©-©©>■ *@a-®e> 


(e) 


B 5 (4) 


o=0 


a=0 


- ViAv uAjJ uAjJ vlAjJ 



FIG. 4: Some basis states in the five sets from B 1 (N P = 4) to 
B 5 {N P = 4). 

^0 + ^elliptic- Then we will treat Vhop in degenerate per- 
turbation theory using eigenstates of Tio + Vcrnptic as t> as i s 
states. 

In both cases we need to find the the eigenstates of 
Ho- We consider only the states with the energies of the 
lowest order of U and comment on the number of the 
degeneracies, g, of each eigenenergy. For N p — 1 the 
eigenenergy and the degeneracy are (E^°\g) = (0,8). 
These basis states B 1 (N p = 1) are shown in Fig. ^ For 
N p = 2 we find (E^°\g) = (V, 16). These basis states 
B 1 (N p = 2) are shown in Fig. For N p = 3 there are 
three groups of basis states B 1 (N P — 3), B 2 (N P = 3) 
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and B 3 (N P = 3). They have {E (0 \g) = {U+2V- J, 24), 
(C7+2V, 16) and (Z7+2F+ J, 8), respectively. For N p = 4 
there are five groups of basis states from B 1 (N p = 4) to 
B 5 (N P = 4). They have (E^,g) = (2U + W - 2 J, 9), 
(2£/ + 4 V - J, 12), (2[/ + 4V, 10), (2U + AV + J, 4), and 
(2£/ + 4F + 2 J, 1), respectively. 

Consider basis states |</>j > that belong to a group 

B l (N p ). Let us call this group I and its energy E^°\ In 
the first order degenerate perturbation theory we diago- 
nalizc the matrix 



(a) 



V; 



(1) 



= < 



\Vr 



pert | 



>, 



(31) 



and find the eigenstates and energies. If this matrix is 
zero we apply the 2nd order degenerate perturbation the- 
ory by diagonalizing the matrix 



y (2) =< 



iv; 



i 



pcrt £(o) 



n 



-v peT t\4>j > 



(32) 



unpert 



Note that the state V per t\(t>j > may not belong to the 
group I; it may belong to another group I' with degener- 
ate energy Ey . The resulting eigenvectors are the linear 
combinations of the basis states \4>j > with appropriate 



expansion coefficients c.y. IV'fe >= X) 



3^1 



bj >. The per- 
turbed eigenenergies are < Vfc 1 Vpert E t _-^ m ert Vpert 1 i>k >■ 
In the parabolic DQDs when the number of electrons is 
odd, the degeneracy of the low energy states for Ho is 
lifted in the 1st order perturbation of Vhop- For a even 
number of electrons, on the other hand, we need to use 
the 2nd order energy. It indicates that the hoppings of 
electrons are more important for the odd-numbered sys- 
tem than for the even-numbered one. All the matrix 
elements Vij have common proportionality factor t or 
t , and consequently the expansion coefficients a in the 
eigenvectors are all independent of these factors. For an 
elliptic potential degenerate perturbation theory is ap- 
plied twice. First we apply it to Vdiiptic, and then apply 
it to Vhop using the eigenstates of Ho + V c iii p tic- 

We have tested our degenerate perturbation theory for 
the following special case of t a = 4 = t, U a = Ub = U 
and U a b = U — J. We have calculated the eigen- 
states, eigenenergies, degeneracies, and quantum num- 
bers. These results are consistent with those for the 
Hamiltonian of the transition metal oxides in Eq. (J2J). 



V. LOW-LYING ENERGY STATES FOR 
PARABOLIC QUANTUM DOTS 

The parabolic DQD is described by the Hamiltonian 
#o + Vhop- Here the low lying energies and the eigen- 
states of our DQD model are investigated by using the 
degenerate perturbation theory for N p < 4 and by trans- 
forming electrons to holes for N p > 5. From now on, in 
order to avoid the long expression, we drop the constant 
term E\N p in the Hamiltonian Ho when writing the en- 
ergy eigenvalues E^> for Ho and E for H. We denote the 



-t 



IVac (1)> 



l¥ac (1)> 



(b) 



i<i (i)>=it 



+ 



FIG. 5: The low-lying energy states for N p = 1 are shown. 
Left and right horizontal bars indicate energy levels of Ho and 
Ho + Vho P - The numbers on the bars indicate the number 
of degeneracy. (a)The energy spectrum and the states are 
shown. Using the particle-hole transformation, we can obtain 
the energy spectrum for N p = 7 from that of N p = 1, see the 
text, (b) One of the groundstates is depicted. 



./Vp-particle eigenfunction for the parabolic Hamiltonian 

(N p ) > where ar and <tt indicate the z- 
component of the total angular momentum and the twice 
of the total spin, respectively. In order to distinguish the 
states with the same values of ax and ctt, we do num- 
bering by assigning k. The state with the lower k value 
has either the same or the lower energy than the higher 
/c-valued one. Since the hopping term can be written as 
Vhop = — 2t£l x we need to diagonalize the first order ma- 
trix elements (Oj,)^ 1 ^. If these first order matrix elements 
are all zero we need to diagonalize the second order ma- 
trix elements (Qx)lj (the definitions of the superscript 1 
and 2 are given in Eqs. I|31|) and H32|) h 

N p = 1: The energy spectrum is depicted in Fig. 0(a). 
Since [Vh op ,^ z ] ^ leads to the broken dot pseudospin 
rotational symmetry, the energy is split into two levels. 
Explicit calculations show that all eigenstates have r = 
1/2, S = 1/2 and = 1/2. Each split level is 4-fold 
degenerate with r z = ±1/2 (a — ±1) and S z — ±1/2 

(a = ±1). For the ground states |Vw > we have fl x — 
1/2. One of the ground states, which has r 2 = 1/2 and 
S z = 1/2, is shown, as an example, in Fig. [31 (b) . Other 
ground states can be obtained by flipping the orbital- 
pseudospin or the spin, or doing both in the example 

state. For the excited states |Vw >, we have il x — 
— 1/2. For all these eigenstates [Ho, t 2 }^^ j0 - T (N p ) >= 
although the operational identity is as [Hq, t 2 ] ^ 0. This 
is because the eigenstates are one electron states. Sec 
Eq. (|16fl and the explanation following it. Consequently, 
t is a good quantum number. For the same reason Q is 
a good quantum number (see Eq. H20|) h The spectrum 
for N p = 7 is exactly the same as that of N p = 1 except 

that the eigenstates corresponding to jf/'io! > have the 
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FIG. 6: The low-lying energy states for N p = 2 are shown: 
(a) the energy spectrum and the states, and (b) two of the 
ground states are depicted. 



energy E = 9U - 3 J + 12V + t while |v4a > have E = 
9U — 3 J + 12V — t. This is because t is replaced by —t 
in the particle-hole transformation. 

N p = 2: In the basis states we include only states with 
one electron in each dot. With these basis states the 
Hamiltonian matrix of Vh op vanishes since the hopping 
term produces states with two electrons on the same 
dot. Therefore we apply the 2nd order degenerate per- 
turbational calculation to find eigenstates. The results 
are shown in Fig. (a). For the ground states, we find 
5=1 and fi = 1. It is 3- fold degenerate with S z = ±1 
and S z = 0. Fig. EJb) shows the graphical representa- 
tions for the ground states with S z — 1 (or = 2) and 
S z = {<jt = 0). The states with S z = — 1 can be ob- 
tained by flipping all spins of the state with S z = 1. Note 
that for all eigenstates each dot contains only one elec- 
tron. Therefore from Eqs. (|16|l and 1)20(1 it follows that 
the numbers r and Q are good quantum numbers of the 
eigenstates. There are competitions between ferromag- 
netic and superexchange interactions in our model Hamil- 
tonian. Two-electron parabolic DQDs are stable when 
the ferromagnetic interactions win, making the spins of 
electrons parallel. On the other hand, the orbital pseu- 
dospins are antiparallcl so that r = 0. 



N„ 



3: Fig. (a) shows how the low energy spec- 



trum for N p = 3 is changed by the perturbation Vhop- 
Note that in the groundstates each dot contains either 
one or two electrons and the dot with two electrons has 
the orbital pseudospin singlet state of Tj = 0. There- 
fore from Eq. I|lfi[) it follows that r is a good quantum 
number of the groundstates. fl also happens to be a good 
quantum number for the groundstates. The ground-state 
energy corresponds to the states with r = |, S = | and 
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±i or S z 



±% and r. 



±4 



\. They have S z - ^ 2 ^ ^ z - ^ 2 , z - 
respectively, yielding the ground state 8-fold degenerate. 
In Fig. [3 (b), we show two ground states of them, one 
with t z — \ and S z = | and the other with t z = \ 
and S z = \. The others can be obtained by flipping all 
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FIG. 7: The low-lying energy states for N p = 3 are shown: 
(a) the energy spectrum and the states, and (b) two of the 
ground states are depicted. 



spins or all orbital pseudospins or both. From the ground 
states, one can notice that our three-electron system fa- 
vors the total spin of electrons to be maximized, showing 
the spin ferromagnetism. On the other hand, the orbital 
and dot pseudospins are minimized. Not all the eigenen- 



(2,3) 



>. 



ergy states are eigenstates of fi , for example, \ipu 
But the numbers 5, S z and r z are good quantum num- 
bers of all the eigenenergy states. 

Np — 4: After the 2nd-order perturbation calculations, 



we obtain the energy spectra in Fig. [S] (a). The bases 
B 1 {A) with the lowest unperturbed energy — 2U — 
2 J + AV are recombined to be the eigenstates with the 
ground-state energy E = 2U-2J + W-6t 2 /(U + J-V), 
and the others with E = 2U - 2 J + W - At 2 /{U + J - V) 

and E = 2U - 2 J + AV . The ground state |^d( 4 ) > is 
unique and has t z = 0, S = and S z = 0. It is also the 
eigenstate for r with r = 0. Note that in the ground- 
states each dot contains two electrons and both dot have 
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by the particle-hole transformation as we mentioned in 
Sec. 11111 The procedure is as below. First, we obtain the 
energy spectra of H h (U, J,V,t; Nh = 8 — N p ), which is 
nothing but the spectra of fi(U, J, V, t; N p = Nh). And 
each eigenstate of TL h is given by simply writing hj aa 
(h!j a(T ) instead of c jacr (c] Qcr ) and |0 >/, instead of |0 > in 
the eigenstate of Ti. Then, we change the hopping energy 
from t to —t in order to obtain the energy spectrum for 
n h (U,J,V,-t;N h ). Finally we add {W - J +W){A-N h ) 
to the energy and transform all states by replac- 
ing h ja(T {h) aa ) with c] QCT (c jaa ) and |0 > h with 

C lll C ll,-l C l,-l,l C l,-l,-l C 211 C 21,-l C 2,-14 C 2,-l,-l|0 > - 

For N p — 8, all of the quantum states on the double 
dots are filled up and the energy of the system is 
E = 12U ~ 4J + 16V". The ground states for N p > 4 
have the same magnetic properties for the states with 
8 — Np. This results from the fact that (2N p — 8) particles 
among N p electrons always remain nonmagnetic, each 
two bound as a pair of electrons with up and down spins 
in the same orbitals on the same dots. For N p = 8, the 
system is naturally nonmagnetic. 



— I <L (4)>. I VZ (4)> 

2U-2]+4V-—tl 

U+J-V 

r 1 d w> 

2U-2J+4V-— °l 

U+J-V 




FIG. 8: The low-lying energy states for N p = 4 are shown: 
(a) the energy spectrum and the states, and (b) the ground 
state is depicted. 



the singlet states of the orbital pseudospins t± = ti = 0. 
Therefore from Eq. (|16l) it follows that t is a good quan- 
tum number of the groundstates. These results suggest 
that the system with N p — 4 is stable when the spins 
and the orbital pseudospins of electrons are antiparallel 
each other. The reason is as follows. For N p = 4, the 
electrons occupy four different quantum states of double 
dots. Unlike the case for N p — 1, 2 and 3, whenever an 
electron of one dot hops onto the other dot it must face 
another electron with the opposite spins due to the Pauli 
exclusion principle. Hence, the superexchange interac- 
tion dominates the ferromagnetic interaction in order to 
lower energy by hopping. It is worth noting that, in the 
ground state, Si = S2 = 1- Within a dot, the spins 
are ferromagnetic. However, Si and S 2 are aligned in 
opposite directions, leading S = 0. Furthermore, the to- 
tal orbital pseudospin r is naturally zero. In the excited 
states, t z , S and S z are good quantum numbers. 

5 < Np < 8: The spectra for N p = 5, 6 and 7 are 
obtained from those for N p = 3, 2 and 1, respectively, 



VI. LOW-LYING ENERGY STATES OF 
ELLIPTIC POTENTIAL 

In the following we will investigate the energy spectra 
for elliptic dots from N p = 1 to N p = 4. We obtain the 
total energies and the eigenstates up to the nonvanishing 
leading terms using two successive degenerate perturba- 
tion calculations, the first, with respect to V^mptic an d 
the second, to Vhop, as we mentioned in Sec. IV. In these 
calculations, the term EiN p in the potential T4mptic as 
well as in the Hamiltonian Tio are ignored. 

We will denote the final eigenstates for 7i to t as 
\^t\n p ) > where ut and k arc cquivalcntly defined 
as in Sec. Note that ar is dropped this time be- 
cause V^niptic breaks the orbital pseudospin's up-down 
symmetry so that the z-component of the total angular 
momentum is not conserved anymore. 

In the first order perturbational calculation with re- 
spect to Kniptic the matrix (V^iiiptic)^- oc (Vj,)-*- has to 
be diagonalized (see Eq. When all matrix elements 

of (V^niptic)^ vanish, the second order matrix elements 

(Kmptic) £ j- is calculated (see Eq. (221) ■ 

N p = 1: The resulting spectra for N p = 1 are shown 
in Fig. El (a) . For N p = 1, the Hamiltonian Ti 

tot is ex- 
actly diagonalized. Unlike the Hamiltonian for parabolic 
dots Tio, the total Hamiltonian does not commute with 
the z-component of the total orbital pseudospin, i.e., 
[Htot, T z ] 7^ 0. This is due to the elliptic potential 14m p tic. 
Hence, after the perturbation Kniptic is introduced, the 
energy is split into two levels. Each level is split again 
into two levels by the hopping term Vh op , which is re- 
sponsible for the fact that TL to t does not commute with 
VL Z , i.e., [7itot,^z] 7^ as we mentioned in the previous 
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FIG. 9: (a)The low-lying energy spectrum for N p = 1 is de- 
picted for an elliptic DQD. The left, middle, and right hor- 
izontal bars indicate energy levels of Ho, Ho + V c m p ti c , and 
Htot = Ho + Kiiiptic + Vkop- (b)One of the groundstates for 
Htot is depicted. 
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FIG. 11: (a)The low-lying energy spectrum for JV P = 3 of an 
elliptic DQD. (b)Two ground states are depicted. 



FIG. 10: (a)The low-lying energy spectrum for N p = 2 of an 
elliptic DQD. (b)The groundstate is depicted. 



section The ground state is depicted in Fig. (b). 
The other ground state can be obtained by flipping the 
electron's spin. For these ground states {tp'S 1 ^ >, we have 
(t,t v ) = (±,-±), (S,S Z ) = (i, \a) and (Q,Q X ) = (1,1). 
Note that in all eigenstates each dot contains either one 
or zero electrons. Therefore from Eqs. (|l(jfl and (|2L)fl it 
follows that r and fl are good quantum numbers. For 
the excited states, all eigenstates have t = S = ^ and 
fi = i. Each energy is doubly degenerate with S z = ±=. 



N„ 



2: The energy spectra of Htot for N p = 2 are 



shown in Fig. (a). Note that similar to the case of 
two-electron parabolic DQDs, the numbers r and il are 
good quantum numbers. The ground state 1^0^(2) > 
is unique and is drawn in Fig. (b). It has {t,t v ) = 
(1,-1), (S = 0,S g ) = (0,0) and Q = 1. One can no- 



tice that the eigenvalues of the total orbital pseudospin 
and the total spin are different from those of the ground 
state of a parabolic DQD with r = and 5=1. Excited 
states |</4 2) (2) > and l^o (3) ( 2 ) > 01 elliptic DQDs have 
the same form as the degenerate groundstates \ipg 1 l (7 (2) > 

and |^q 1 q(2) > of parabolic DQDs. In these ground 
states of parabolic DQDs, the total orbital pseudospin of 
two electrons forms a singlet state, (| 1UI-> — | JJ.f|->)/\/2 
with t = 0. Since Vciiiptic flips all pseudospins in turn, 
this singlet state vanishes as 2eT y ^(\ lH!-> — | JJ-1r>) = 

-ie-^fl U> +\ U> -I n> -I U>) = 0. Hence, the 
antisymmetric property of the orbital pseudospin singlet 
leads to the unchanged total energy. In the groundstate 
of a two-electron elliptic DQD, the total spin is mini- 
mized while the total orbital pseudospin is maximized, 
showing orbital pseudospin ferromagnetism. 

N p = 3: The results of the elliptic DQD for N p = 3 
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/2U+4V+2J+- 



l<*'(4)> 



TABLE I: Groundstate magnetic properties of parabolic 
DQDs for different values of N p . A hyphen means that the 
relevant quantum number is not a good quantum number. 



are shown in Fig. 1111 (a). Each unperturbed energy is 
split into two levels by V^mpti C : the lower energy state 
has t v = — g and the higher energy state t v = +■=. 
They are split again when the hopping perturbation 
Vh p is introduced, as shown in the spectra in the right 
side. The ground states are four-fold degenerate with 
S z = ±i and S z = ±|. Two of them are depicted 
in Fig. [H] (b), and others can be obtained by flip- 
ping all spins in the depicted states. In the ground 
states, we have (t 7 t v ) = (|,— |), S = | and (fl,£l x ) = 
{\,\)- For the same reason as in the parabolic DQD 
for N p = 3, r and are good quantum numbers in 
the ground states. Similar to the parabolic DQD, the 
ground states show the spin ferromagnetism while the 
orbital and the dot pseudospins are minimized. This can 
be understood as follows: the ground states for ellip- 
tic dots are the linear combinations of the parabolic-dot 
ground states with the same spin quantum number <tt, 
i.e., |^(3) >= ^(1^1(3) > -#^,3,(3) >) and 

|^ (1) (3) >= ^(|^(3) > -#^(3) >). Hence, it is 
natural that they show the same magnetic properties of 
the ground states of parabolic dots. Note that the num- 
bers S and S z are good quantum numbers in the excited 
states. 

N p = 4: Fig. E| (a) shows the energy spectra for 

N p = 4. The energy of = 2U + 4V-2J and the cor- 
responding basis states B 1 (4) are not changed after the 
perturbation Vdiiptic- The reason is as follows. Each dot 
in the basis states has Sj = 1, j = 1,2 (triplet states). 
The spins of two electrons within a dot are parallel which 
does not permit the change of the orbital pseudospins of 
the electrons due to Pauli exclusion principle. Hence, 
the elliptic potential does not change the basis states 
B 1 (4). Since only the hopping term Vh op has an effect on 
those states, the ground state is the same as that for the 
parabolic dots. We depict it in Fig. EH(b). Of course, 
5 = and r = in the ground state. In the excited 
states, S and S z are, in general, good quantum numbers. 

5 < Np < 8: The elliptic potential is unchanged under 
the particle-hole transformation. The energy spectra and 
the states for N p > 5 can thus be obtained from those 
for N p < 4 via the particle-hole transformation. 
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VII. CONCLUSION AND DISCUSSION 
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FIG. 12: (a)The low-lying energy spectrum for N p = 4 of an 
elliptic DQD. (b)The ground state is depicted. 



We have investigated theoretically low excitation states 
in the 'p' molecular orbital of parabolic and elliptic 
DQDs. In our many body calculation we include intra- 
and inter-dot electron-electron interactions, in addition 
to the intradot exchange interaction of 'p' orbitals. Wave- 
functions, energy splittings, and degeneracies of states 
of the molecular orbitals are determined. The low ly- 
ing excited states may be probed experimentally by 
measurin g cu rrent when a finite source drain voltage 
is present^, ^| . Our result shows how many current 
peaks should be observed in such an experiment. 

We have also determined the properties of orbital mag- 
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TABLE II: Groundstate magnetic properties of elliptic DQDs 
for different values of N p . 

netism of groundstates. Since the magnetization is di- 
rectly related to the z-component of the total angular 
momentum our results may be used to determine the 
value of the magnetization of DQDs[12j. Note that the 
z-component of the total angular momentum is in fact 
equal to the value of the z-component of the total pseu- 
dospin. In elliptic DQDs the z-componcnt of the total 
angular momentum is not a good quantum number. 

The orbital pseudospin and spin properties for 
parabolic and elliptic DQDs are summarized in Tables 
I and II, respectively. We see for parabolic DQDs that 
the total spin of electrons takes the maximum possible 
value except for N p = 4 and N p = 8. For these values 
of N p the total spin is minimum. The orbital pseudospin 
is zero for N p = 2,4,6,8. Note that for both N p = 2 



and 6 the spins and pseudospins are ferromagnetic and 
antiferromagnetic, respectively, i.e., S — 1 and r = 0. 
However, for elliptic DQDs the opposite is true. In this 
case the groundstate energy is minimized when the total 
spin takes the value S — 0. This difference in the in- 
terplay between orbital pseudospin and spin magnetism 
in parabolic and elliptic DQDs can be tested experimen- 
tally. When a weak magnetic field is applied the ground- 
state energy of the parabolic DQD should split into three 
while that of the elliptic dot remains unsplit. A sufficient 
deformation from the parabolicity is required since the 
strength of the ellipticity assumed to be stronger than 
the magnitude of the tunneling energy in our results. In 
our calculations we have assumed that the exchange in- 
teraction between dots is negligible. In some cases the 
interdot exchange interaction can be important, and can 
lead to canted phases j^- Further investigations includ- 
ing the interdot exchange are needed to test the stability 
of canted states in our DQDs. 
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